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Materials which can undergo extremely fast displacive transformations as well as very slow 
diffusive transformations are studied using a Ginzburg-Landau framework to understand the physics 
behind microstructure formation and time-temperature-transformation (TTT) diagrams. This 
simple model captures the essential features of alloys such as steels and predicts the formation 
of mixed microstructures by an interplay between diffusive and displacive mechanisms. The 
intrinsic volume changes associated with the transformations stabilize mixed microstructures such 
as martensite-retained austenite (responsible for the existence of a martensite finish temperature) 
and martensite-pearlite. 
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Based on kinetic considerations, solid-to-solid phase 
transformations can be broadly categorized as diffusive 
and displacive Diffusive phase transformations - 
such as spinodal decomposition — involve long-range mo- 
tion of atoms and are therefore slow. In displacive phase 
transformations — e.g. martensitic transformations |2| — 
the crystal structure changes through a unit cell distor- 
tion, without any long-range motion of the atoms. 

These two types of phase transformations may interact 
or compete, e.g. in steel 0, Ti-Al-Nb H, Cu-Al-Ag Q, 
Pu-Ga , and NiTi shape memory alloys Q • Eutectoid 
steels decompose into pearlite, i.e. ferrite plus cementite 
(iron carbide, Fe3C), by a diffusive process. However, 
at low temperature carbon diffusion is very slow and, 
although it is not the ground state, martensite can form 
on fast cooling. At intermediate temperature a different 
microstructure, bainite, may form in steel through a 
mixed diffusive-displacive mechanism j|. Information 
on these phase transformations and their kinetics can 
be conveniently represented on time-transformation- 
temperature (TTT) diagrams. A theoretical understand- 
ing of the interplay between diffusive and displacive 
modes is paramount to comprehend TTT diagrams and 
microstructures, which are important technologically in 
the design of heat treatments. 

Olson et al. |£j theoretically studied coupled diffusive 
and displacive transformations to calculate a TTT dia- 
gram for bainite. However, this one-dimensional model 
cannot describe the microstructural complexity and the 
long-range elastic fields generated during such transfor- 
mations. The Ginzburg-Landau method (also known as 
the phase-field method) on the other hand can naturally 
describe the appropriate habit planes without assuming 
any given microstructure. It has been extensively used to 
study microstructural evolutions in diffusive E3 El EM 
and martensitic 0, 0, EH EH E3 transformations, 
but (apart from one work on a TTT diagram for ferrite 
and martensite in iron EH) these phase transformations 
have typically been studied separately. 



Using the Ginzburg-Landau theory, we study a model 
system which can undergo both a phase separation and 
a square-to-rectangle martensitic transformation. We 
label "pearlite" any phase-separated region, even if the 
microstructure does not correspond to real pearlite per 
se. Although this letter focuses on TTT diagrams, 
other thermal treatments are possible using this model, 
e.g. continuous cooling transformation (CCT). Likewise 
martensitic transformations other than the square-to- 
rectangle studied here can be modeled. 

The free energy of the system is expressed as 
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Here g e \ is the usual non-linear elastic free energ y de nsity 
for a square-to-rectangle martensitic transition [lj, [l^ : 
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where e\ = (s xx + e y y)/V2 is the hydrostatic strain, e 2 = 
(?xx — £yy)/V% is the deviatoric strain, and es = e xy 
is the shear strain. The {£ij} are the linearized strain 
tensor components. T is the temperature and Tm and Tp 
are constants pertaining to the austenite-martensite and 
austenite-pearlite phase transformations respectively, c 
is the composition (c = corresponds to the composition 
of austenite and the composition of pearlite is c = =Lcq). 
The chemical free energy is given by E3 
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and a coupling between elastic distortions and composi- 
tion is introduced as g cp \ = X2 C c 2 e2- Notice that in the 
stress- free conditions, at equilibrium e\ = x\ c c + x±2 e\\ 
during the transformations, volume changes may be 
introduced, either due to a lattice mismatch between 
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FIG. 1: The bulk energy as a function of e2 and c at (a) T — 
0.2, (b) T = 0.7, and (c) T — 1.5. Darker areas correspond 
to lower energies. 



the phase-separating components {x\ c ) or due to the 
transformation strains (#12). 

Since we are interested in a qualitative understanding 
of the physical mechanisms, we choose simple values for 
the parameters: A2 = 2, A4 = 4, Aq = 9.6, B\ = 1, 
£3 = 1, C 2 = 6, C 4 = 12, x 2c = 5, T M = 0.5, 
Tp = 1, fc c = 2, and /c e 2 = 0.1. The coupling constants 
#12 and xi c are varied in different cases to understand 
the effect of volume changes. The homogeneous part 
of the free energy is depicted in Fig. as a function 
of 62 and c at different temperatures. The different 
phases can be identified as follows: austenite corresponds 
to (c = 0, e2 = 0), martensite to (c = 0, e2 7^ 0), and 
pearlite to (c 7^ 0, e2 = 0). Above Tp only austenite is 
stable, Fig. |l(c)| Between Tm and Tp austenite and 
martensite are unstable and pearlite is the ground state, 
Fig. |l(b)| Below Tm, pearlite is the ground state and 
martensite is metastable, Fig. |l(a)| 

The kinetics of the transformations are described by 
equations of motion for the displacement fields and the 
composition. The evolution of the displacements is 
described by [HE} 
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where S is the functional derivative, p is a density, and 
v is the time derivative of the displacements, u. The 
second term on the right-hand side in Eq. (|3a|) is a viscous 
damping term. The strains, {sij}, are then obtained 
as derivatives of the displacements. The evolution of 
the composition is described through the Cahn-Hillard 
equation 

dc(r,t) , « SG 
at oc(r, t) 

where M is the temperature-dependent mobility, 

M = M exp(-Q/T). (5) 

Note that this temperature dependence of this mobility 
is necessary in order to obtain realistic TTT diagrams. 
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FIG. 2: (color) (a) TTT diagram (A: austenite, M: marten- 
site, and P: pearlite) and the microstructures at (b) T — 0.9 
and t = 2 000, (c) T = 0.5 and t = 2 000, and (d) T = 0.49 
and t = 50 for X12 — x\ c — 0. Red and yellow: martensite; 
light and dark blue: pearlite. 



All simulations are two-dimensional with periodic 
boundaries and on a 128 x 128 lattice (5x = 1 and 
St 0.2). We use rj 0.01, p = 1, M = 2, and 
Q = 5. The initial system, made of 100% austenite, 
includes random fluctuations around c = and u = 0. 
The system is quenched instantaneously to temperature 
T and held at this temperature. For each value of T 
we record the times at which 10% martensite and 10% 
pearlite form, as well as the time at which the austenite 
content drops below 10%. 

Figure [2(a)] shows the resulting TTT diagram for the 
X12 = %ic — case. The austenite-pearlite phase 
transformation requires diffusion and therefore time. At 
low temperature, diffusion is slow and so is pearlite 
formation. Martensite on the other hand forms through 
a displacive mechanism, which does not require long- 
range motion of atoms. Figure |2(a)| shows that at 
temperatures below about T = 0.5 (the "martensite 
start temperature") austenite transforms to martensite 
(Fig. |2(dj| ) and pearlite forms only at higher temperature. 
At temperatures close to Tp = 1 the undercooling is 
so small that pearlite formation is very slow. Pearlite 
formation is also delayed at temperatures close to T = 
0.5 because of the slow diffusion. Consequently there 
exists an intermediate temperature at which pearlite 
formation is the fastest and a C-curve can be observed 
around T = 0.8 in Fig. |2(a)| The typical length scale 
of pearlite depends on diffusion length and therefore 
on temperature: at low temperatures pearlite is fine 
— Fig. |2(c)| — and it is coarser close to Tp — Fig. |2(b)| — , 
consistent with experimental observations. Note that the 
amounts of the two components of pearlite are equal 
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FIG. 3: (color) Microstructures for X12 = and x\ c = at 
T = 0.49 and (a) t = 4 000, (b) t = 6 000, and (c) t = 13 000. 
Red and yellow: martensite; light and dark blue: pearlite. 



because their compositions are symmetric with respect 
to that of austenite (Fig.[Q. 

Along with Fig. |2(d)[ Fig. |3| shows the microstructure 
evolution at T = 0.49 (just below the martensite start 
temperature). At t = 50 only martensite is found, 
Fig. |2(d)| At t = 4 000, pearlite has already started 
forming at the interfaces between the martensite variants 
aligned along [1 1] and those along [1 1] , Fig. |3(a)[ that 
is at interfaces which are relatively high in energy. For 
longer times pearlite grows at the expense of martensite 
(Fig. |3(bJ| ) and the latter finally disappears (Fig. |3(c)| ). 
This appearance of pearlite in the martensite region 
occurs at progressively later times at lower temperatures. 

As pearlite is the ground state and martensite is 
only metastable the latter forms only because pearlite 
formation is slow at this temperature. It is therefore 
expected that pearlite can eventually form in a marten- 
sitic system as this will decrease the energy. However 
there exists an energy barrier because, unlike austenite, 
martensite is metastable. This makes pearlite nucleation 
from martensite slower than nucleation from austenite 
and accounts for the discontinuity of the 10% pearlite line 
across T = 0.5 in Fig. |2(a)] Note that the model cannot 
be expected to accurately estimate the time necessary for 
pearlite to form as the nucleation process is not properly 
described. 

So far we assumed that there was no bulk strain 
associated with pearlite or martensite formation, i.e. 
#12 = %ic = in Eq. (JJ). We now relax this constraint 
to study the effect of the term x±2 e\ of Eq. JIJ (keeping 
xi c set to zero). This term couples deviatoric strain (i.e. 
martensite) and hydrostatic strain: it associates a net 
volume change with the martensitic transformation. The 
effect of X12 on the microstructure at T = 0.49 can be 
seen in Fig. El Whereas for low values of X12 austenite 
completely transforms to martensite (e.g. Fig. |4(a)| ), 
for larger values of X12 martensite formation results in 
stresses so large that a system made of pure martensite 
would be unstable. Consequently, part of the system 
remains austenitic, Fig. 14(c)] A similar microstructure 
was obtained by Onuki [16]. This retained austenite 
transforms to pearlite, which then grows at the expense 
of martensite (Fig. |4(cjJ| ) . 

The term x\ c c in Eq. (jT]). similar to Ref. 10. corre- 
sponds to a lattice mismatch between the two compo- 




FIG. 4: (color) Microstructures at T = 0.49. Top row: x\ c = 
0, bottom row: x\ c — 1; two left columns: x\2 — 0, two right 
columns: x\2 = 1; columns 1 and 3: t = 50, columns 2 and 
4: t = 50 000. Red and yellow: martensite; green: austenite; 
light and dark blue: pearlite. 



nents of pearlite (i.e. ferrite and cementite in the case of 
steel), such that the average lattice spacing of pearlite is 
that of austenite. Figure^lshows the effect of X12 and x\ c 
on the microstructure. Since x\2 can stabilize austenite, 
at short times its value determines the microstructure: 
Figs. |4(a)| and |4(e)| (which correspond to xyi — 0) 
show pure martensite whereas in Figs. |4(c)| and |4(g)| 
(x\2 — 1) there is a mixture of martensite and austenite. 
The long-time microstructures on the other hand are 
controlled by the value of x\ c \ in Figs. |4(b)| and |4(d)| 
(corresponding to x\ c = 0) there is only pearlite whereas 
Figs. |4(f)| and |4(h)| [x\ c = 1) show a mixed pearlite- 
martensite microstructure (which does not disappear in 
longer simulations). x\ c stabilizes martensite at long 
times: when x\ c = 1 a purely pearlitic system would 
not be stable because of the large stresses it would 
generate. This is a kind of stress-induced martensitic 
transformation. 

As expected |2l| martensite-martensite interfaces are 
oriented along (11). The orientations of martensite- 
austenite and martensite-pearlite interfaces are different. 
For an interface at an angle 6 with respect to [10], 
elastic compatibility requires that cos 20 = 5e\/5e2 where 
Se± and Se2 are the variations of e\ and e2 across the 
interface. The values of e\ and e2 in Fig. |4(g)| give 
w ±30° or « ±60°, which is consistent with 
the interface orientations observed in the figure. In 
Fig. |4(h)| the rounded shape of pearlite grains aims at 
minimizing interface energy; therefore the orientation of 
the interfaces is not determined by elasticity alone. 

Figure focuses on x\2 = %ic = 1- Several features 
of the TTT diagram shown in Fig. |5(a)| are different 
from the TTT diagram obtained for x±2 = x\ c = 0, 
Fig. |2(a)| There exists a martensite finish temperature 
(around T = 0.45) due to retained austenite. Martensite 
can be found above T = 0.5 because — due to underlying 
dilation strains and elastic compatibility — a mixture 
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FIG. 5: (color) (a): TTT diagram for xi 2 = xi c = 1 (A: 
austenite, M: martensite, and P: pearlite). Microstructures 
at T = 0.49: (b) t = 1 000, (c) t = 1 500, and (d) t = 2 000 
(red and yellow: martensite; light and dark blue: pearlite). 



of pearlite and martensite is more stable than pure 
pearlite (this remains true up to T w 0.73). This 
splitting of the pearlite region in the TTT diagram is 
akin to the pearlite-bainite transition in steel (however, 
the microstructure shown in Fig. |4(h)| is different from 
that of bainite). When x±2 = x\ c = 0, below T = 0.5 
pearlite can nucleate heterogeneously at the interface 
between martensite variants; this process is rather slow. 
When x\2 = x\ c = 1, even below T = 0.5 some 
austenite remains and pearlite formation proceeds in this 
retained austenite. Since austenite is unstable there is 
no nucleation barrier and there is no discontinuity of 
the 10% pearlite line (the pearlite nucleation is slower 
at lower T only because the diffusivity is lower). 

Figures |4(g)[ |5(b)H5(d)| and |4(h)| show the time 
evolution of the microstructure at T = 0.49 for x\2 — 
x\ c = 1. By t = 1 000, pearlite nucleated in the retained 
austenite and new martensite formed where there was 
none before (pearlite cannot form alone because of 
the stresses it generates), Fig. |5(b)| After austenite 
disappears pearlite grows at the expense of the large 
grains of "primary" martensite and makes them "rotate" 
towards (01), Fig. |5(c)1 From large martensite grains 
alternating with areas completely devoid of martensite 
at t = 50, Fig. |4(g)| the system evolves to a state where 
martensite is more homogeneously distributed, Fig. |5(d)| 
at t = 2 000, through a double mechanism of martensite 
formation and martensite destruction, and finally to a 
coarsened microstructure, Fig. |4(h)| The feature size 
thus goes from large to small to medium with increasing 
time. Consequently the "final" structure is independent 



of the initial one, both in terms of grain size and interface 
orientation (compare Fig. |4(h)| to Fig. |4(g)| ). 

A new Ginzburg-Landau approach has been proposed 
to study alloys — such as steels — which can undergo dis- 
placive as well as diffusive transformations. It captures 
the important features of TTT diagrams and microstruc- 
tures and sheds some light on the role of the interplay 
between the two types of transformations in stabilizing 
mixed microstructures. The existence of a martensite 
finish temperature (i.e. of retained austenite) is due to 
an hydrostatic strain associated with the martensitic 
transformation. When a strain is associated with pearlite 
formation, martensite and pearlite coexist at intermedi- 
ate temperatures, i.e. in the region of the TTT diagram 
where bainite is typically found in steels. The model 
also shows that in these mixed microstructures the habit 
planes can be different from the pure martensite case. 

We thank David Srolovitz and Srikanth Vedantam for 
useful discussions. 
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